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Abstract 

A planetary microlensing event occurs when a planet perturbs one of the two 

images created in a point-mass microlensing event, causing a deviation from the 

standard Paczyhski curve. Determination of the two physical parameters that can 

be extracted from a planetary microlensing event, the planet/star mass ratio q, 

and the planet/star separation in units of the stellar Einstein ring, y p , is hampered 

by several types of degeneracies. There are two distinct and qualitatively differ- 
• <— i , 

^ . ent classes of planetary events: major and minor image perturbations. For major 

image perturbations, there is a potentially crippling continuous degeneracy in q 
which is of order 87 , where 5^ is the maximum fractional deviation of the plan- 
etary perturbation. Since the threshold of detection is expected to be 5^ ~ 5%, 
this degeneracy in q can be a factor of ~ 20. For minor image perturbations, the 
continuous degeneracy in q is considerably less severe, and is typically less than a 
factor 4. We show that these degeneracies can be resolved by observations from 
dedicated telescopes on several continents together with optical/infrared photom- 
etry from one of these sites. There also exists a class of discrete degeneracies. 
These are typically easy to resolve given good temporal coverage of the planetary 
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event. Unambiguous interpretation of planetary microlensing events requires the 
resolution of both types of degeneracy. We describe the degeneracies in detail 
and specify the situations in which they are problematic. We also describe how 
individual planet masses and physical projected separations can be measured. 

Subject Headings: gravitational lensing - planetary systems 

submitted to The Astrophysical Journal: October 15, 1996 
Preprint: OSU-TA-26/96 
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1. Introduction 



Two world-wide networks are currently searching for extra-solar planetary sys- 
tems by making densely sampled observations of ongoing microlensing events to- 
ward the Galactic bulge (PLANET, Albrow et al. 1996; GMAN, Pratt et al. 1996). 
Several other groups will join the search shortly and there is serious discussion of 
new initiatives that would intensify the search by an order of magnitude. More than 
100 microlensing events have been detected to date by four groups, MACHO (Al- 
cock et al. 1996), EROS (Ansari et al. 1996), OGLE (Udalski et al. 1994), and DUO 
(Alard 1996) based on observations made once or twice per night. The events typ- 
ically last one week to a few months. MACHO and OGLE have reported "alerts", 
events detected before peak. This alert capability is what has allowed PLANET 
and GMAN to make intensive, sometimes round-the-clock, follow-up observations 
in hopes of finding the planetary perturbations which are expected to last a day 
or less. 

In sharp contrast to this explosion of observational activity, theoretical work 
on planet detection has been rather sparse, amounting to only five papers in as 
many years. Mao & Paczyhski (1991) originally suggested that planets might 
be detected in microlensing events. Gould & Loeb (1992) developed a formalism 
for understanding the character of planetary perturbations and made systematic 
estimates of the rate of detection for various planetary-system parameters. Bolatto 
& Falco (1994) studied the detection rate in the more general context of binary 
systems. These early works assumed that the lensed star could be treated as a point 
source. The usefulness of this approximation depends primarily on the angular size 
of the source 9*, relative to the planetary Einstein ring, 9 P , 



Here 6 e is the Einstein ring of the lensing star, m and M are the masses of the planet 
and its parent star, and D Q \, D\ a , and D os are the distances between the observer, 
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lens, and source. For Jupiter-mass planets at typical distances (D\ s ~ 2 kpc) from 
bulge giant sources, 8 p ~ 30* so the approximation is a reasonable one. However, 
for Saturns, Neptunes, and especially Earths, the finite size of the source becomes 
quite important, and even for Jupiters it is not completely negligible. Moreover, as 
we will stress below, it is quite possible to mistake a "Jupiter event" in which the 
source size is negligible for a "Neptune event" with > 9 P . Hence it is essential 
to understand finite-source effects even to interpret events where the source size is 
in fact small. 

Progress on finite-source effects was substantially delayed by problems of com- 
putation. Like all binary lenses, planetary-systems have caustics, curves in the 
source plane where a point-source is infinitely magnified as two images either ap- 
pear or disappear. If one attempts to integrate the magnification of a finite source 
that crosses a caustic, one is plagued with numerical instabilities near the caus- 
tic. While it is straight forward to solve these problems for any given geometry, 
the broad range of possible geometries makes it difficult to develop an algorithm 
sufficiently robust for a statistical study of lensing events. Bennett & Rhie (1996) 
solved this problem by integrating in the image plane (where the variation of the 
magnification is smooth) rather than the source plane (where it is discontinuous). 
They were thereby able to investigate for the first time the detectability of Earth 
to Neptune mass planets. Gould & Gaucherel (1996) showed that this approach 
could be simplified from a two-dimensional integral over the image of the source to 
a one-dimensional integral over its boundary. The implementation of this method 
requires some care. We describe the practical procedures elsewhere (Gaudi 1996). 
The difficult computational problems originally posed by finite-source effects are 
now completely solved. 

To date, the analysis of planetary-system lensing events has focused on the 
question of "detectability" which was quantified by Gould & Loeb (1992) as a 
certain minimal fractional deviation from a standard Paczyhski (1986) light curve 
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where x is the projected lens-source separation in units of 6 e . Note that this curve 
is characterized by just three parameters: to the time of closest approach, /3, the 
impact parameter in units of 9 e , and t e , the Einstein radius crossing time. Bennett 
& Rhie (1996) adopted a similar approach but added the qualification that the 
deviation persist for a certain minimum time. 

Here we investigate a different question: How well can the parameters of the 
planetary-system be measured? As discussed by Gould & Loeb (1992), if there 
are light-curve data of "sufficient quality", two planetary-system parameters can 
generically be extracted from a microlensing event that displays a planetary per- 
turbation. These are the planet/star mass ratio, q, and the planet/star projected 
separation in units of the stellar Einstein ring, y p , 



Here a p is the physical projected separation, and r e = D o \0 e . As we discuss in 
§ 8, it will often be possible to make additional observations that specify the 
mass and distance of the lensing star, or equivalently M and r e . For these cases, 
the measurements of q and r e yield the mass m = qM and projected separation 

a p = Vp r e- 

If a planet were detected by observing a deviation from the standard curve, but 
its mass ratio remained uncertain by a factor of 10, the scientific value of the de- 
tection would be severely degraded. Indeed such "detections" would probably not 
receive general acceptance. Thus, the problems of planet detection and parameter 
measurement are intimately connected. Microlensing planet-detection programs 
must monitor a total of at least several hundred events in order to obtain repre- 
sentative statistics on the frequency of planets. These observations require large 
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blocks of 1-2 meter class telescope time coordinated over several continents. For 
funding agencies and time allocation committees to make rational decisions about 
the allocation of scarce resources, and for observers to make rational choices among 
prospective targets, it is essential to determine what are the minimum observational 
requirement for detecting planetary systems and measuring the characteristics of 
the detected systems. 



2. Types of Degeneracy 

2.1. Discrete 

Planetary-system lensing events are subject to two different discrete degenera- 
cies. The first ambiguity relates to which image the planet is perturbing: the major 
image outside the Einstein ring or the minor image inside the Einstein ring. For 
almost all cases, this degeneracy is easily broken provided there is good temporal 
coverage of the light curve. However, if it is not broken the uncertainty in q and 
y p can be a factor of a few. The magnitudes of these uncertainties depend only on 
the overall geometry of the event and not on the mass of the planet. The second 
ambiguity relates to whether the planet lies closer to or farther from the star than 
does the position of the source image that it is perturbing. This degeneracy is 
more difficult to break, but it does not seriously affect the determination of q, and 
the uncertainty induced in y p is proportional to g 1 / 2 and is therefore often much 
smaller than the one induced by first degeneracy. These two discrete degeneracies 
are illustrated in Figure 1. The values of q and y p for each of the four possible 
solutions are displayed in Table 1. 
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Figure 1. Discrete degeneracies. Panel (a) shows a lensing light curve with (solid curve) and without (dashed 
curve) taking account of the presence of a planet with mass ratio q = 10 -3 . Panel (b) shows the associated lensing 
geometry. The two solid curves represent the path of the images relative to the lens. The crosses represent the 
image positions at the time of the perturbation. The circles are the four planet positions for which the light curves 
reproduce the measured parameters 5d (maximum fraction deviation) and t& (FWHM of deviation) at the peak of 
the disturbance when the source-lens separation is Xd- The filled circle is the "actual" planet position. Panel (c) 
shows the four associated light curves for times near the peak of the perturbation, i ,<z- Note that time is expressed in 
units of the perturbation time scale, td, not t e . The bold curve corresponds to the "actual" planet position. Clearly, 
if the light curve is well sampled, the two dashed curves corresponding to the image position inside the Einstein ring 
in panel (b) could be ruled out immediately. However, the two solid curves are less easily distinguished. These differ 
by ~ 15% in planet/star separation and 10% in mass. 3ee panel (b) and Table 1. 



TABLE 1 



Degenerate Parameter Values: Discrete 



planet/star planet/star 

separation mass ratio 

Vp q/% 

Major Image 1.40 1.00 

1.19 0.91 

Minor Image 0.75 1.08 

0.80 0.88 



2.2. Continuous 

In addition, there is a continuous degeneracy arising from finite-source effects 
being misinterpreted as a larger value of q. This is because q is determined from the 
(square of the) duration of the planetary perturbation relative to the total duration 
of the event. If the size of the source is larger than the Einstein ring of the planet, 
then the duration of the planetary perturbation will be the crossing time of the 
source, not of the planet Einstein ring. Figure 2 shows 10 light curves all with 
the same maximum fractional deviation, 5d, and same full width half maximum 
(FWHM) of perturbation, td- The parameter that differs in each of these curves is 
the ratio of source radius, to planet Einstein radius, 9 P = q9 e , 



Table 2 gives the inferred values of q and of the proper motion fi (of the planetary 
system relative to the observer-source line of sight) associated with each curve in 
units of the "fiducial" values associated with p = 0.3. In so far as one could not 
distinguish among these curves, any of these parameter combinations would be 
acceptable. The fiducial parameters qo and /xq would then be measurable but the 
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Figure 2. Ten light curves with shear 7 = 0.6, <j> = 90° (see eqs. 3.1 and 3.2), all with maximum 
deviation Sd = 10% and FWHM td — 0.06 t e . The ratios of source radius to planet Einstein ring 
range from p = 0.1 to p = 2.87, the largest source radius consistent with this maximum deviation. 
Table 2 gives the corresponding values of q = m/M, and proper motion, p, relative to the fiducial 
values q and p at the arbitrarily chosen value p = 0.3. 



actual values of \i and q would not. The proper motion of both bulge and disk lenses 
is typically ji ~ 

0(^lsr/^o) ~ SOkms^kpc- 1 , where Vlsr ~ 220 km s" 1 is the 
rotation speed of the Local Standard of Rest, and Rq ~ 8 kpc is the Galactocentric 
distance. If, for the example shown in Table 2, the fiducial value were measured as 
Ho ~ ^lsr/^O; one might then choose to argue that the proper motions associated 
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with the low-mass solutions (i.e. /i ~ Ato/3) would be so low as to be a priori 
unlikely. However, these solutions could not actually be ruled out by such an 
argument, since the distribution of fi is rather broad (see Han & Gould 1995). 
Thus, there would remain a factor ~ 15 uncertainty in the planet/star mass ratio. 

TABLE 2 



Degenerate Parameter Values: Continuous Major Image 



U-llilcllolUlllCoo 


JJldllL,l / oldl 






llldob IdllO 




P 


q/qo 


fl/flO 


0.10 


1.095 


2.867 


0.20 


1.041 


1.470 


0.30 


1.000 


1.000 


0.60 


0.957 


0.511 


0.90 


0.767 


0.381 


1.20 


0.566 


0.332 


1.50 


0.373 


0.327 


1.80 


0.236 


0.343 


2.10 


0.163 


0.354 


2.40 


0.127 


0.351 


2.70 


0.093 


0.364 


2.87 


0.074 


0.383 



2.3. Relation Between Degeneracies in q and /i 

From the relation /i = 9 e /t e , we obtain the identity = (0 e /9p)(9 p /6*)(6*/t e ) 

or 

/W /2 = r- (2.2) 

Since the quantities on the right hand side of this equation are observables, the 
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product on the left hand side must be constant for all allowed parameter com- 
binations in any given planetary event: ppq 1 / 2 =constant. This equation then 
establishes a relationship between degeneracies in q and degeneracies in p. If a 
range of solutions are permitted that have different values of p but very similar 
values of q, then we say that the mass ratio is not degenerate. However, it follows 
from equation (2.2) that the proper motion p then varies inversely as p and there- 
fore that it is degenerate. Similarly, if the range of allowed solutions all have the 
same value of p, then the proper motion is not degenerate but then q oc p~ 2 and 
so the mass ratio is degenerate. This relationship is illustrated by Table 2. The 
region p < 0.3 has well-determined q but degenerate p, while the region p <; 1 has 
well-determined p but degenerate q. 



3. The Chang-Refsdal Lens Approximation 

In order to systematically investigate the role of these degeneracies and to 
determine the data that are required to break them, we follow Gould & Loeb (1992) 
and approximate the planetary perturbation as a Chang-Refsdal lens (Chang & 
Refsdal 1979; Schneider, Ehlers, & Falco 1992). A Chang-Refsdal lens is a point 
mass (in this case the planet) superimposed on a uniform background shear 7. 
For any given lensing event, the value of 7 is simply the shear due to the lensing 
star at the unperturbed position of the image that is perturbed by the planet. The 
evaluation of 7 is made at the mid-point of the perturbation. The source position 
at this midpoint, x d is known from the light curve (see Figs. l-a,b). The associated 
image positions are yd ± with shears t± given by 



2 (4 + A) 1 / 2 ±x d 
7± = V d ,±, Vd,± = 7, • (3-1) 



Thus the shear is known (up to a two-fold ambiguity) simply from the position of 
the planetary perturbation on the overall light curve. 
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When computing the fractional deviation, 5, of the Chang- Refsdal lens from the 
standard Paczynski curve, one always normalizes relative to the total unperturbed 
magnification [eq. (1.2)] which includes both the image perturbed by the planet 
and the image that remains unperturbed (Gould & Loeb 1992). 

The Chang-Refsdal approximation permits an immense conceptual simplifica- 
tion of the problem. For a point source, all possible light curves of an event with a 
given Xd can be represented on a pair of diagrams, one for 7+ and one for 7_. All 
possible planetary perturbations can therefore be represented by a single-parameter 
family of such diagrams. See Figure 3. For a given event, one knows j3 and Xd 
from the overall light curve. One can therefore compute 7± using equation (3.1) 
and thereby pick out which two diagrams are relevant. One also knows the angle 
(f> at which the source cuts through the diagram, 

sin0 = — . (3.2) 

If, for example, Xd = 0.516 and (3 = 0.4, then all possible light curves are repre- 
sented by the parallel lines indicated in Figures 3-c,d. If the light curve is well 
sampled, it is easy to distinguish between 7 + and 7_. Suppose that 7+ (Fig. 3-b) 
is correct, and say that the maximum fractional deviation is 5d = 20%. Then one 
can immediately identify the correct curve as being either B or D. The observed 
duration of the perturbation relative to that of the whole event then sets the scale 
of diagram relative to 9 e and thus determines the mass ratio. 

Of course, one does not know a priori that finite source effects can be ig- 
nored. However, for any Xd, all possible events can still be represented by a 
single-parameter family of diagrams. The relevant parameter is, p, the ratio of 
the angular radius of the source to the Einstein radius of the planet. Hence, it is 
quite easy to study all possible degeneracies. See Figure 4. 

The drawback of using the Chang-Refsdal approximation is that it is not exact. 
Moreover, for any given lensing event, it is straight forward to construct models 
that are exact. As we argue below, however, the lack of exactness has no significant 
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Figure 3. Chang-Refsdal magnification contours of a point source as a function of source position in units of the 
planet Einstein ring, 9 P = q6 e , for various pairs of shears (7+ and 7_ = 7+ 1 ) corresponding to planetary perturbations 
of the major and minor images, respectively. Magnification contours are calculated including the contribution of the 
unperturbed image. Contour pairs are for 7+ = 0.8, 0.6, 0.4, and 0.2 in panels (a,b), (c,d), (e,f), and (g,h), and 
correspond to source positions at perturbation of x d = 0.22, 0.52, 0.95, and 1.79. Super-bold contour is no deviation. 
Bold contours are S = 5%, 10%, 20%, and 00. Non-bold contours are —5%, —10%, and —20%. Diagonal lines 
in panels (c) and (d) represent possible trajectories assuming that the overall light curve shows Xd = 0.52 (i.e. 
7+ = 0.6) and (3 = 0.4 (i.e. </> = sin -1 fi/xd ~ 50°). If the maximum deviation were observed to be 5d = 20% (and 
the point-source approximation were known to be valid), then the trajectory must be either B or D. 
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7=0.60 



7=1.67 




Figure 4. Chang-Refsdal magnification contours of a non-point source as a function of source position in units of 
Op, for the pair of shears (7+, 7-) = (0.6, 1.67) corresponding to planetary perturbations at source position Xd — 0.52. 
The ratios of planet Einstein radius to source radius are p — 0.5, 1.0, 1.5, and 2.5. Compare with point-source case 
shown in Figs. 3-c,d. Contours levels are the same as for Fig. 3. 
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impact in the analysis of degeneracies. On the other hand, using the exact solution 
increases the dimensionality of parameter space and thereby the conceptual com- 
plexity of the problem, without any compensating benefits. We therefore strongly 
advocate using the Chang-Refsdal framework. We present a more detailed analysis 
in the Appendix. 

4. Degeneracy Between Major and Minor Images 

Figure 3 indicates that it should generally be quite easy to distinguish between 
perturbations of the major and minor images provided that there is good tempo- 
ral coverage: perturbations of the major image have one major positive excursion, 
while perturbations of the minor image have two positive excursions separated by a 
large negative excursion. However, if the observations are made from only one site, 
then good temporal coverage is far from automatic. The time scale for these excur- 
sions is the minimum of the crossing time of the star, 9*//i ~ 10 hours for a giant, 
and the crossing time of the planet Einstein ring 9 p / jj, ~ 10 (m/50M©) 1//2 hours. 
Thus it would be quite possible to observe a positive excursion (or a significant 
fraction of it) one night and then miss any subsequent excursions due to one or 
two nights of bad weather. However, if there were three observing sites on different 
continents, such large data gaps would be rare. 

There is nevertheless another possible source of degeneracy between the major 
and minor images. If <fi is sufficiently small, then a source coming close to one of 
the caustics of the perturbation of the minor image could cross the star-planet axis 
at a point far enough from the planet that the negative excursion along this axis 
would be very small. In this case, the light curve might be mistaken for one due 
to a perturbation of the major image. In Figure 5, we present examples of this 
degeneracy for three different values of 7+ for = 0, along with the corresponding 
color shifts A(V — H) (see § 6.3). The parameter p was chosen in each case 
such that the curves for 7+ and 7_ would be most similar. It is clear that the 
degenerate curves could be distinguished only if precise measurements could be 
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Figure 5. Major/minor image degeneracy for (j> = 0, with 7+ = 0.61,0.38,0.25, in panels a,c,and e. 
solid curve corresponds to 7+, the dashed curve to 7_ = 7+ ■ Also shown in panels b,d, and f are 
associated fractional color change A(V — H). 



made at the wings of the perturbation. Considering now the curves for A(V — H), 
there are relatively large (2 — 10%) fractional color changes associated with the 
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7+ curves throughout the event, while the fractional color changes associated with 
the 7_ curves are always negligible. This large difference in the magnitude of the 
fractional color change between the 7+ and 7_ curves arises from the fact that the 
gradient of the magnification across the face of the star is much larger in the 7+ 
case, and thus the color effects are therefore more pronounced (see § 6.3). Thus 
by measuring A(V — H), one can distinguish between the two degenerate cases 
even during the peak of the perturbation. The larger the value of 0, the larger the 
negative excursion in the minor-image perturbation, and thus curves with > 
will be less degenerate than the examples shown in Figure 5. 



5. Degeneracy of Planet Position Relative to Unperturbed Image 

In general, this degeneracy introduces an uncertainty in y p which is Ay p ~ 
2a9 p / 9 e where a6 p is the separation between the planet and the unperturbed im- 
age at the mid-point of the perturbation. From Figure 3 one sees that if per- 
turbations 5d ~ 5% are detectable, the typical planetary event will have a ~ 5. 
Hence, if the degeneracy remains unbroken, the fractional uncertainty is Ay p /y p ~ 
(lOOm/M) 1 / 2 . For Jupiter -mass planets in orbit around M dwarfs, this error is of 
order unity, while for Neptune-mass planets it is ~ 10%. On the other hand, the 
degeneracy in q is small (see Table 1 and also Appendix). 

In this section we consider only point sources. If the planet has a low mass, so 
that finite source effects are important, then (as mentioned above) the difference 
between the two degenerate solutions is small and distinguishing between them 
is relatively less important. In addition, finite-source effects are more properly 
addressed in the context of the mass/finite-source degeneracy discussed in § 6. 
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5.1. Perturbations of the Major Image 

It is clear from Figure 3 that it is impossible to break the degeneracy if = 
90°, i.e., if the planetary perturbation takes place at the peak of the light curve 
(xd = (3). If the source crosses the perturbation structure in the region of the 
caustic (a ^ 1) then degeneracy is relatively unimportant and, in any event, is 
easily broken (provided < 90°) due to the richness of the structure in this region. 
We therefore focus on the case a > 1. From Figure 3, one sees that there is an 
asymmetry in the light curve which has opposite senses depending on whether the 
source crosses to the left or the right of the planet. If it passes to the right, then 
the deviation is more pronounced at the beginning of the perturbation than at the 
end, and if it passes to the left, the deviation is more pronounced at the end. We 
define the asymmetry factor P^ as the maximum over all times t of the fractional 
difference, 

P ( 7 , S d ) = max{|5(t , d + t) - 6(t , d - t)\} (5.1) 

where to,d is the mid-point of the perturbation and S(t) is the fractional deviation 
as a function of time. To lowest order, one may approximate 

P ( 7 ,5 d )=P(7,^)cot0. (5.2) 

From Figure 3, one can see that as 7 increases, the positive contours of 5 become 
more stretched along the planet-star axis, and thus low-peak perturbations occur 
farther from the areas of negative excursion. One would therefore expect smaller 
values of P for larger values of 7. Figure 6 shows P as a function of 5d for several 
values of 7. As expected, P generally decreases with increasing 7. From Figure 6, 
one can determine, for a given sensitivity, whether the degeneracy can be broken 
for any perturbation. For example, consider trajectories such that <fi ^ 75°. If 
one were sensitive to asymmetries of P^ ~ 1%, i.e., P ~ .04, then the degeneracy 
could be broken only for 7 ^ 0.3 , and then only if 5^ > 0.1. On the other hand, if 
~ 30°(P ~ 0.004), then the degeneracy could be broken for essentially all values 
of 7. 
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Figure 6. Asymmetry factor P (see eqs. (5.1) and (5.2)) for nine values of 7, as a function of maximum 
deviation, 64. The actual asymmetry of a given light curve, P<j>{^(, Sd), is given by P^ij, Sd) ~ -P(7, Sd) cot 4>. 
Using this figure and formula, one can therefore determine whether the degeneracy can be broken for any 
given sensitivity threshold. 

5.2. Perturbations of the Minor Image 

As we discussed in § 4, to distinguish minor-image from major-image pertur- 
bations, it is necessary to observe the negative excursion (centered on the of 
the right-hand side of Fig. 3. If these are observed, then one can easily distinguish 
the case where the source transits the right side from the case where it transits 
the left side of the x-axis (provided < 90°). Hence, there is no degeneracy of 
minor-image perturbations unless the more severe major /minor-image degeneracy 
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remains unbroken. 



6. Continuous Mass Degeneracy of Major Image Perturbations 

By far the potentially most crippling form of parameter degeneracy is the one 
that is illustrated in Figure 2 and is tabulated in Table 2. The basic character 
of this degeneracy can be understood analytically using the following theorem 
(Gould & Gaucherel 1996): if the unperturbed major image crosses the position of 
the planet and the source is larger than the major-image caustic structure, then 

s ^7Wr Ah) = r^S- (6 - 1) 

The FWHM of such an event is td ~ 2(csc 0)pg 1 / 2 t e . On the other hand, the 
FWHM of a low-peak perturbation of a point-source is td ~ 2(csc <p)q 1 / 2 t e . Suppose 
that an event has observables t e , j3, Xd, td, and 5d- One can form the combination 
of observables Q = [{(3 / Xd){td/t e ) /2] 2 , and can obtain one possible solution that 
reproduces the maximum deviation and FWHM: 



q~Q, p<l. (6.2) 



However, the solution 



Q (2 1-7 2 \ 1/2 , , 

q ~7u: "~ pm - a "" = Utt?J • (63) 

would also reproduce the height and width of the curve. Note that the ratio 
of masses for the two solutions is p 2 nax . For 5d ~ 5%, this ratio is typically <; 
20 and can be as high as 40. Thus, unless this degeneracy is broken, any low- 
peak perturbation of a point source by a Jupiter-mass planet can masquerade as a 
Neptune- mass event, and vice versa. All intervening masses are permitted as well. 
Clearly, unless this degeneracy is broken, low-peak perturbations will contain very 
little information about mass, and unambiguous detection of low-mass planets will 
be impossible. There are three possible paths to breaking this degeneracy. 
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Figure 7. Fractional deviations 6 for i7-band light curves2Bimilar to Fig. 2, except now shown for shears 7 = 0.2, 0.4, 0.6, 
and 0.8, and for trajectory of source motion <f> = 30°, 45°, 60°, and 90° (left to right). Corresponding curves for <p —* —<t> can 
be found by reversing the x-axes. In each curve, the maximum deviation is 84 = 10% and the FWHM is td — 0.06 t e . 



6.1. Proper Motion Measurement 

If the proper motion p of the lens relative to the source were measured, then 
one could partially break the degeneracy between equations (6.2) and (6.3). The 
angular size of the source 9* would be known from its dereddened color and mag- 
nitude and Stefan's law. The time for the source to cross the perturbation x-axis, 
t c ~ 29* csc (f)/p, would then also be known. If one found t c < tj, this would im- 
ply that td was dominated by the size of the planet Einstein ring, not the source. 
Hence, the solution (6.2) would be indicated. On the other hand, if t c ~ td, one 
would know only that the solution (6.2) was not correct, and that the mass ratio 
lay somewhere in the interval p m a X <3 < 1 < Q- See Table 2. 

Proper motions can be measured if the lensing star transits or nearly transits 
the source or by imaging the split image of the source using infrared interferometry. 
See Gould (1996) for a review. Another approach would be simply to wait a few 
decades and measure the angular separation of the lens and source. Since typical 
proper motions are 5 mas yr _1 , the separation should be ~ 0.1" after a few decades. 
Unfortunately, most lenses are probably fainter than Mj = 10, while typical giant 
sources are Mj ~ 0, and even turnoff stars are Mj ~ 3. Thus it would be difficult 
to image the lens until it was quite well separated from the source at which point it 
might be hard to distinguish it from random field stars. We explore this possibility 
further in § 8. 

6.2. Detailed Light Curves 

Although the parameter combinations (6.2) and (6.3) reproduce the gross fea- 
tures of the perturbation (peak and FWHM) equally well, the detailed structures 
of the light curves are different. Figure 2 illustrates the principal difference for 
elongated perturbation structures, in this case 7 = 0.6. When p ~ p max , the wings 
show a dip because the source passes over the caustic which is surrounded by re- 
gions of negative perturbation (see Figs. 3-c and 4). On the other hand, when 
p ^ 1 the approach to the peak is smooth because the source is passing over the 
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smooth outer portion of the ridge seen in Figure 3-c. From Figure 2 it is clear that 
if these wing structures could be resolved at the ~ 1% level, then the degeneracy 
in mass could be reduced from the factor ~ 15 seen in Table 2, to a factor ~ 1.5. 

Figure 7 is an array of 16 diagrams each similar to Figure 2, but with different 
values of 7 (0.2, 0.4, 0.6, and 0.8) and different angles of source motion (30°, 45°, 
60°, and 90°). It is clear that it is easier to break the degeneracy as 7 increases 
and as <fi decreases. For 7 ^ 0.4 the uncertainty in q could be significantly reduced 
if the wing structures could be resolved at the ~ 1% level. For 7 ^ 0.2, however, 
distinguishing between the degenerate curves would require an accuracy <C 1%. 
Figure 8 shows the values of q/qo and /i//io as a function of p for each of the 
combinations of and 7 in Figure 7. Note that larger values of p m ax are allowed 
for smaller values of 7, and thus the range of acceptable values of q is largest for 
small 7. This is especially disturbing in light of the fact that the degenerate curves 
are most similar for small 7. 

6.3. Optical/Infrared Colors 

A major shortcoming of the detailed-light-curve method for breaking the de- 
generacy is that it depends critically on obtaining accurate observations during 
two brief intervals covering the wings of the light curve. As a practical matter, it 
may be difficult to obtain such coverage for a variety of reasons. Once the event 
is noticed, observatories that are dedicated to the planet search can engage in 
frequent monitoring and thereby obtain very accurate light curves. However, it is 
quite possible, indeed likely, that the planetary perturbation will not be recognized 
in time for intensive monitoring of the first wing. Sometimes observation of the 
first wing is crucial to breaking the degeneracy. Moreover, the second wing will 
likely be observable from at most one observatory which could be affected by bad 
weather. 

Optical/infrared color measurements by contrast yield degeneracy-breaking in- 
formation throughout the event. The reason is that by the principle of equivalence, 
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Figure 9. Fractional color change A(V — H) for light curves shown in Fig. 7. 



lensing of point source is achromatic. If lensing introduces color changes, the lens 
must be resolving the source (Witt 1995; Loeb & Sasselov 1995). The best op- 
portunity to observe this effect is by looking for optical/infrared color differences 
(Gould & Welch 1996) because giant stars are more limb-darkened in the optical 
than in the infrared (Manduca, Bell & Gustafsson 1977). Thus, if the planet Ein- 
stein ring is larger than the source (and the low peak is due to the source passing 
over regions of small perturbation), the color changes will be very small. On the 
other hand, if 6* > 9 P (and the low peak occurs when the large source passes over 
the caustic), the caustic structure will resolve the differential limb-darkening of the 
star and the color changes will be more pronounced. Figure 9 shows the V — H 
colors for the same parameters as are used for the H band curves as in Figure 
7. The magnitude of the fractional color change is largest for smallest 7. This 
is fortunate, since, as discussed in § 6.2, the degeneracy is most severe for small 
7, both in terms of the similarity in the light curves, and in the range of allowed 
values of q. It is therefore essential to have optical/infrared color measurements to 
ensure that the continuous degeneracy can be broken for all possible values of 7. 

7. Continuous Degeneracy of Minor Image Perturbations 

There is also a continuous degeneracy for minor-image perturbations, but the 
degeneracy is considerably less severe than for major-image perturbations because 
the caustic structure is qualitatively different. As with major image perturbations, 
the basic character of the minor image degeneracy can be understood analytically. 
Consider the following theorem (Gould & Gaucherel 1996): if the unperturbed 
minor image crosses the position of the planet and the source encloses both minor- 
image caustics, then 

where /(p, 7) = [(1/2 + p~ 2 ) 2 - 7 2 p~ 4 ]~ 1//2 , and the limit applies for p > 7 1 / 2 . 
That is, in contrast to the major-image perturbation [cf. eq. (6.1)], the minor-image 
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perturbation goes to zero rapidly for large sources. For minor image perturbations, 
the caustics are located at (Schneider et al. 1992) 

4aus ~ 2( 7 - I) 1 ' 2 . (7.2) 

Thus the source must have p ^ 2(7 — l) 1 / 2 to enclose both caustics. If the source 
is significantly larger than this, the perturbation will be negligibly small. Hence 
we can restrict attention to sources p < d cans . 

For minor-image perturbations, is the FWHM of the negative deviation. For 
a low-peak point source perturbation (a ^ 2), this duration scales as the distance 
between the contours of 5 = at a. (See Fig. 3). For finite sources of increasing 
p, the trajectories must move closer to the center to maintain the observed value 
of 5d- (See Fig. 4). However, until p becomes so large as to cover the caustics, the 
positions of the 5 = contours basically do not change. (Compare Fig. 3d with Fig. 
4). Since these contours are approximately horizontal, td is not greatly influenced 
by changes in p. Finally, since the largest permitted source has p ~ <i caus , which 
is also approximately equal to the separation of the 5 = contours, td remains 
roughly the same even for this extreme case (see Fig. 3). The small degeneracy 
that does exist arises from the difference between this extreme case on the one 
hand and the smaller sources and point sources on the other. Examining Figure 
3, we can expect that the degeneracy will be somewhat larger for larger values of 
7, since the contours of 5 = become less horizontal as 7 increases. In Table 3 we 
give the degeneracy in the inferred values of q for 5d = —10% and —5%, and for 
several values of 7. Note that the largest degeneracy in q is only a factor of ~ 4. 
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TABLE 3 



Continuous Minor Image Degeneracy 



mass ratio 
degeneracy 



5 d 


7 


Qmax./ Qmin 


-10% 


1.25 


1.45 




1.43 


1.09 




1.67 


1.68 




2.00 


3.82 


-5% 


1.25 


1.66 




1.43 


1.16 




1.67 


1.22 




2.00 


1.67 




2.50 


2.34 



Figure 10a shows twelve light curves for 7 = 1.67 and = 90°, all with 
maximum negative perturbation 5d = —10%, and all with the same FWHM. Table 
4 gives the inferred values of q and p for each curve, relative to the fiducial values 
go and po associated with p = 0.3. Note that the degeneracy in the derived mass 
ratios is only a factor of ~ 1.5. Also note that the inferred mass ratios of the 
first nine curves agree to ~ 4%. Thus to resolve the small degeneracy in q, one 
only needs to distinguish between the last four curves. From Figure 10a it is clear 
that this would be possible if one could resolve the positive perturbation structures 
at the ~ 1% level. Furthermore, the situation presented in Figure 10, for which 
4> = 90°, is the worst case scenario. Due to the structure of the caustics of minor- 
image perturbations, trajectories with < 90° display marked asymmetry about 
to,d, excepting trajectories with a ~ 0, which are nearly symmetric. This enables 
one to distinguish between curves with a ^ 1 and a ~ more easily when < 90°. 
This is demonstrated in Figure 11, which shows twelve light curves with the same 
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Figure 10. a) Ten light curves with 7_ = 1.67, <j> — 90°, all with maximum deviation 5d — —10% and 
FWHM td — 0.06 t e . The ratios of source radius to planet Einstein ring range from p = 0.1 to p = 1.83, 
the largest source radius consistent with this maximum deviation. The corresponding relative values of 
q = m/M, and the relative proper motion, p,/p,o, ar e given in Table 4. Source radii of p = 1.6, 1.7, 1.83, and 
1.80 correspond to the bold, dashed, bold dotted, and bold dashed curves, b) The fractional color change 
A(V — H) for the ten curves in panel (a). 29 
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Figure 11. Ten light curves with 7- = 0.6 -1 , <fi = 60°, all with maximum deviation Sd = —10% and 
FWHM td = 0.06 t e . All other parameters are the same as Fig. 10, and are given in Table 4. Curves are as 
in Fig. 10. 

parameters as Figure 10, except that now = 60°. Comparing Figs. 10 and 11, it is 
clear that the curves are appreciably less degenerate for = 60° than for <fi = 90°. 
From Figure 10b, we see that the magnitude of the fractional color change for 
perturbations with 7 = 1.67 is always small, A(V — H) ^ 1%. From numerical 
calculations, we find that A(V — H) & \% regardless of the value of 7. Thus, 
in contrast to major image perturbations, optical/infrared color are not useful in 
resolving the degeneracy in minor image perturbations, since the magnitude of the 
fractional color change is always small. 
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TABLE 4 



Degenerate Parameter Values: Continuous Minor Image 



impact 


dimensionless 


planet / star 


proper motion 


parameter 


source radius 


mass ratio 




OL 


o 

r 


a /an 


111 Lift 


4 32 


1 


QQ3 


4 014 


4 30 


20 


1 016 

X . W X W 


1 984 


4.25 


40 


1 000 

X .www 


1 000 

X .WWW 


4 21 


60 

w . w w 


1 006 

X .www 


66 s ) 

w .WWW 


4.00 


0.80 


1.020 


0.495 


3.70 


1.00 


1.033 


0.394 


3.50 


1.20 


1.017 


0.331 


3.20 


1.40 


0.994 


0.287 


2.60 


1.60 


1.004 


0.249 


2.00 


1.70 


1.093 


0.225 


1.00 


1.83 


1.591 


0.173 


0.00 


1.80 


1.628 


0.174 



8. From Mass Ratios to Planet Masses 

If the various degeneracies described in this paper are broken, one generally 
recovers two planetary-system parameters from a planetary micro lensing event: q 
and y p . While q is of some interest in its own right, y p is not. The quantities one 
would most like to know are the planet mass m = qM and the physical projected 
separation a p = r e y p . One could take a purely statistical approach to estimating 
these quantities: given the measured time scale t e of the event and a plausible 
model of the distribution and velocities of lenses and sources along the line of 
sight, r e and M can be estimated to a factor of 3. In this section, we discuss what 
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further constraints might be obtained on M and r e in order to determine m and 
a p . 

The single most powerful method of acquiring additional information would be 
to launch a parallax satellite (Refsdal 1966; Gould 1995a; Gaudi & Gould 1997) 
which would routinely measure f e = (D os / D\ s )r e and often measure the direction 
of motion as well. This information would, by itself, narrow the uncertainty in the 
mass to a factor ~ 1.7 (see Han & Gould 1995, especially figure 7). However, if 
the proper motion fx were also measured, this would yield a complete solution of 
the lensing geometry including both M and r e (e.g. Gould 1996). In general, one 
expects to measure fx only in ~ 20% of giant events even with relatively aggressive 
observations (Gould 1996). However, for events with planetary perturbations, fx 
can be measured much more frequently. Recall from § 6 that for giant sources, 
major image perturbations, and planetary masses m ^ 100 M®, the planet usually 
resolves the source (if it is detected at all) and that in the process of resolving 
the resulting degeneracy, one measures fx. Even when < 9 p , the source will 
sometimes cross a caustic in which case fx can be measured. Finally, for 6* < p 
one can obtain a lower limit fx > fx m i n based on the lack of detection of finite- 
source effects. Since the mass is given by M = (c 2 /AG)f e t e fx, and since f e and t e 
are measured, this gives a lower limit on the mass (Gould 1995b). 

However, it is much more difficult to resolve the finite source degeneracy for 
minor-image perturbations. Even though (or rather, because) the measurement of 
the mass ratio, q, is not seriously hampered by this degeneracy, the proper motion 
fx is poorly determined. See § 2.3. Thus, it would be necessary to measure fx using 
other methods. See § 6.1. 

We now address several questions related to one of those methods: direct imag- 
ing of the source and lensing star several decades after the event. For definiteness, 
we suppose that the measurement is made after 20 years. The expected separation 
is ~ 0."1, but could plausibly be ~ O.'^. At Baade's Window, the expected number 
of stars Mj < 10 inside this radius is ~ 0.5 (Light et al. 1996). Thus one would not 
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be overwhelmed with candidates. On the other hand, the great majority of lensing 
events are almost certainly due to objects that are fainter than Mj = 10 simply 
because one does not come close to accounting for the observed events from the 
observed (Mj < 10) stars alone (Han 1996). Thus, to positively identify a candi- 
date star as the lens, one needs additional information. A parallax satellite could 
provide two pieces of corroborating data. First, the measured r e together with the 
proper motion inferred from the candidate-source angular separation would give a 
mass and distance to the lens (Gould 1995b). One could then predict an apparent 
magnitude and see if it agreed with that of the candidate. Second, if the parallax 
measurement gave the angle of motion, one could check this against the direction 
of the source-candidate separation vector. In addition, the candidate's inferred 
proper motion must satisfy the lower limit derived from lack of finite-source ef- 
fects as discussed above. Finally, one could wait another decade or so to see if the 
direction of the candidate's proper motion was indeed away from the source. 
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APPENDIX Justification for the Chang- Refsdal Approximation 

What errors are introduced by the Chang-Refsdal approximation? The unper- 
turbed image structure consists of two images separated by > 29 e . The planet, 
with an effective sphere of influence ~ 9 P <C 6 e can have a major effect on it at most 
one of these. For definiteness, say this is the major image. In the Chang-Refsdal 
approximation, the minor image is then treated as being completely unaffected by 
the presence of the planet. In fact, the planet will change the shear at the minor 
image by 0(6 p /6 e ) and therefore change the magnification by a similar amount. 
However, what is directly of interest for analyzing the planetary perturbation is not 
the absolute difference in magnification with and without the planet. Rather, it is 
the change in this difference over the lifetime of the planetary perturbation. Hence, 
the net effect is 0[(9 p /9 e ) 2 ], i.e., of higher order than the effects being analyzed. 
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We now turn to the errors in the Chang- Refsdal estimate of the magnification 
of the perturbed image. In general, the perturbed image is split by the planet 
into two or four images. For each such image, i, the shear due to the parent 
star is 7$. If this value were exactly equal to 7, the shear at the position of 
the unperturbed image, then the Chang- Refsdal approximation would be exact. 
Typically, A7; = 7^ — 7 is small, A^i/^y ~ 0(9 p /9 e ), so one expects that the errors 
induced by the approximation are small. 

We focus first on perturbations of the major image. Let A9 be the separation 
between the planet and the unperturbed image and define a = A9/9 P . Consider 
first the case a ^> 1 which is important when 7 <: 0.5 because the magnification 
contours then become significantly elongated (see Fig. 3). The image is then split 
into two images, one very close to planet and the other very close to the unper- 
turbed image. For the image close to the planet, the shear due to the parent star 
may be significantly misestimated, A^i/j ~ a9 p /9 e . However, for this image, the 
total shear is dominated by the planet and is 0(a 2 ), so true fractional error is only 
~ a~ 1 9 p /9 e . Moreover, the magnification of this image is small, (9(a~ 4 ), so the 
total error induced by the approximation is ~ a~ 5 9 p /9 e and is completely negli- 
gible. The other image is displaced by ~ oT x 9 p from the unperturbed image, so 
A'ji/'j ~ oT x 9 p j9 e which induces a similar small change in magnification. Recall 
from Figure 3-c, that the source trajectory is determined up to a two-fold degen- 
eracy from the maximum magnification. Since the sign of the image displacement 
is different for the two allowed solutions, the error in estimating the magnification 
structure could result in two types of errors. First, there is an error in the planet 
star separation, but this is only ~ aT x 9 p and is therefore lower by oT x than the ba- 
sic degeneracy indicated in Figure 3-c. Second, there is an error in the estimate of 
q and, in fact a degeneracy because the error has opposite sign for the two allowed 
solutions. This could in principle be significant because, within the Chang-Refsdal 
framework, the two allowed solutions indicated in Figure 3-c have identical values 
of q, and this effect is therefore the lowest order degeneracy. However, the mass 
ratio is estimated from the FWHM of the light curve which is only a weak function 
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of position along the elongated magnification contours. Moreover, the misestimate 
of that position is small. We therefore estimate a fractional mass degeneracy of 
Aq/q ~ a~ 2 q 1 / 2 . 

For a ^ 1 and sources that are small compared to the caustic structure (seen 
e.g., in Fig. 3), the situation is similar to that of caustic-crossing binary- lens events. 
The light curves are highly non-degenerate, and one determines not only q and x p , 
but also p. From the standpoint of understanding degeneracies, the important case 
is when the source is of order or larger than the caustic. Here, there are roughly 
equally magnified images displace roughly by 9 P on either side of the planet. Hence, 
the lowest order errors cancel and the next order errors are ~ 0[(9 p /9 e ) 2 }, and can 
therefore be ignored. 

There is one exception to this conclusion. In the argument given above, we 
implicitly assumed that the planetary perturbation would be significant only over 
an interval of source motion ~ 9 p . This assumption fails when the perturbation 
structure is elongated (7 ^ 0.5) and when the angle of source motion is low (sin <fi = 
(3/xd <C 1). In this case, the local shear is no longer well approximated by the 
shear at the center of the perturbation. A proper calculation would then require 
that the shear be recalculated at every point along the source trajectory, holding 
the planet fixed. This was the approach of Gould & Loeb (1992) and the resulting 
magnification for fixed planet position can be seen in their Figure 3. [In the present 
work, by contrast, what is held fixed in constructing Figs. 3 and 4 is the observable: 
the shear at mid-point of the perturbation.] As can be seen by comparing Figure 
3-c of Gould & Loeb (1992) and Figure 3-c of the present work, for Jupiter mass 
planets the difference in contours can be significant. However, there are three 
points to note. First, such events are rare both because the conditions (7 <; 
0.5, (3 <C Xd) together imply (3 ^ 0.2 and because the elongated contours are 
encountered "edge on", so the cross section is only ~ 9 p /9 e . Second, the effect is 
proportional to q 1 ! 2 and so would not be significant for, e.g., Earth- mass planets. 
Third, the nature of the effect is to provide information to break degeneracies in 
cases when the Chang-Refsdal approximation would lead one to believe that there 
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is no information. In brief, in certain rare cases, the Chang-Refsdal approximation 
leads one to underestimate the amount of information available. 

For perturbations of the minor image, the two principle sources of degeneracy 
are first, confusion of the two caustic peaks with each other and second, confusion 
of one of these peaks with a perturbation of the major image. Because these peaks 
are offset in the direction perpendicular to the star-planet axis, the error in their 
location is 0[(9 p /9 e ) 2 } and hence of higher order than their separation. As in the 
case of the major image, there are certain rare events with sin 1 for which the 
Chang-Refsdal approximation makes the degeneracy seem somewhat worse than it 
is. 



36 



REFERENCES 

1. Alard, C. 1996, in Proc. IAU Symp. 173, Astrophysical Applications of 
Gravitational Lensing, p. 214, (Eds. C. S. Kochanek, J. N. Hewitt), Kluwer 
Academic Publishers 

2. Albrow M., et al. 1996, in Proc. IAU Symp. 173, Astrophysical Applications 
of Gravitational Lensing, p. 227 (Eds. C. S. Kochanek, J. N. Hewitt), Kluwer 
Academic Publishers 

3. Alcock, C, et al. 1996, ApJ, submitted 

4. Ansari, et al. 1996, A&A, in press 

5. Bennett, D., & Rhie, H. 1996, ApJ, in press 

6. Bolatto, & Falco, E. 1994, ApJ, 436, 112 

7. Chang, K., & Refsdal, S. 1979, Nature, 282, 561 

8. Gaudi, B. 1996, in preparation 

9. Gaudi, B., & Gould, A. 1997, ApJ, 477, 000 

10. Gould, A. 1995a, ApJ, 441, L21 

11. Gould, A. 1995b, ApJ, 447, 491 

12. Gould, A. 1996, PASP, 108, 465 

13. Gould, A., & Gaucherel, 1996, ApJ, 477, 000 

14. Gould, A., & Loeb, A. 1992, ApJ, 396, 104 

15. Gould, A., & Welch, D. 1996, ApJ, 464, 212 

16. Loeb, A., & Sasselov, D. 1995, ApJ, 449, 33L 

17. Manduca, A., Bell, R. A., & Gustafsson, B. 1977, A&A, 61,809 

18. Mao, S., & Paczyhski, B. 1991, ApJ, 374, 37 

19. Paczyhski, B. 1986, ApJ, 304, 1 

37 



20. Pratt et al. 1996, in Proc. IAU Symp. 173, Astrophysical Applications of 
Gravitational Lensing, p. 221 (Eds. C. S. Kochanek, J. N. Hewitt), Kluwer 
Academic Publishers 

21. Refsdal, S. 1964, MNRAS, 128, 295 

22. Schneider, P., Ehlers, J., & Falco, E. E. 1992, Gravitational Lenses (Berlin: 
Springer- Verlag) 

23. Udalski, A., et al. 1994, Acta Astronomica, 44, 165 

24. Witt, H. 1995, ApJ, 449, 42 



38 



